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fH ' Abstract. One class of compact binaries of special interest is that of double white 

dwarfs (DWDs). For many of these systems, the exact nature of the evolutionary chan- 
I ■ nels by which they form remains uncertain. The canonical explanation calls for the pro- 

genitor binary system to undergo two subsequent mass-transfer events, both of which 
are unstable and lead to a common envelope (CE) phase. However, it has been shown 
that if both CE events obey the standard oqe prescription, it is not possible to repro- 
duce all of the observed systems. As an alternative prescription, the y-formalism was 
proposed, which parametrizes the fraction of angular momentum carried away with 
mass loss, in contrast to the ckce prescription, which parametrizes energy loss. We 
demonstrate that the y-prescription is also inadequate in describing the evolution of 
an arbitrary DWD binary; clearly we require a deeper understanding of the physical 
mechanisms underlying their formation. We then present a detailed model for the evo- 
lution of Red Giant - Main Sequence binaries during the first episode of mass transfer, 
and demonstrate that their evolution into DWDs need not arise through two phases of 
dynamical mass loss. Instead, the first episode of dramatic mass loss may be stable, 
non-conservative mass transfer. The second phase is then well described by the acE 
prescription. We find that the considered progenitors can reproduce the properties of 
the observed helium DWDs in which the younger component is the more massive. 
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1. Problems with Parametrizing CE Evolution 

In order to explain the dramatic loss of mass and angular momentum needed to form the 
majori ty of observed close binaries with at least one degenerate component, Paczynski 
19761) and lOstrikeri (Q976) first suggested such systems undergo a common envelope 



(CE) phase during their evolution. In this event one component expands on the giant 
branch to engulf the other leading to spiral-in and the removal of the envelope. In the 
standard picture, the consequent orbital shrinkage is parametrized in a simple manner 
in terms of the mass lost and an efficiency factor acE (though see Ivanova, this volume): 
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dLivio & SokedH988r) where m a is the mass of the companion, a x and Of are the initial 
and final separation, r c i is the radius of the donor at the onset of the event, and m&, m^ c , 
and /Md.e are the masses of the donor, its core, and its envelope, respectively. 

However, problems soon became apparent in this straight-forward picture, most 
glaringly in attempting to apply this prescription to the study of double white dwarfs 
(DWDs). Such binaries are considered one of our best means of constraining CE evo- 
lution, since due to their close orbit they must necessarily have passed through a CE as 
their most recent episode of mass loss. At the same time, their small size (the majority 
are found as low-mass helium dwarfs) means they could not possibly have formed by 
any normal process of single-stellar evolution within a Hubble time, implying both 
components must have come about through substantial binar y interaction (i. e . two 
phases of dramatic mass loss while on the first giant branch). Nelemans et all (l2000l) 
demonstrated that, were such systems to have evolved through two CE phases, many 
would have needed to expel the initial donor's envelope with a negative efficiency (»ce) 
in the first event. At the same time, an initial phase of stable, conservative case B mass 
transf er appeared to be an unlikely scenario dNelemans et al.ll2~000l : Ivan der Sluys et al" 



2006). This at the time suggested that an alternative mechanism may be needed in order 
to explain such systems. 

In response to the apparent need for a negative efficiency in the first phase of 
dramatic mass loss, iNelemans et all ( 2000h suggested that rather than guessing at the 
energetics of the event, one could invoke an alternative parametrization in terms of the 
angular momentum lost: 
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where the parameter y describes the average specific angular momentum of the binary 
carried away by the mass shed from the system. This allowed for a significant shrinkage 
of the orbit in the second episode of mass loss with only a modest reduction or even 
wi dening of the orbit in the first phase, dependent on the input masses. This in turn 
led INelemans & Tout! (l2005h to claim that a range 1.5 < y < 1.75 could satisfactorily 
reproduce all post-CE binaries, for specific ranges of initial masses in each case. 

However, while fits for y in this range were found to be common to all those sys- 
tems used in their calibration, there remains quite a wide range in plausible values of 
y for each indi vidual system, corresponding to a range of possible pro genitor masses. 
Prior studies ( Webbink 120081 : Izorotovic et al.ll2O10l : Davis et aT1l2010l also see Zoro- 
tovic et al, and Davis et al, this volume) have demonstrated that in general the use of 
the y formalism rather poorly constrains the possible set of solutions in reconstructing 
post-common-envelope binaries. 

The reasoning behind the great multiplicity of solutions offered by such a narrow 
range of y values becomes clear if we examine the effect of varying y for a fixed system. 
(Note that the discussion here follows that of I Woods" et all (1201 lh . § 3). Consider the 
parametrization given in equation |2l if we substitute for the initial and final angular 
momenta: 
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then we can solve for the initial over the final separation: 
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Figure 1. Effect of small changes in y on the final period, for 2.0 + 0.5M Q (left) 
and 2.0 + 1.5M binaries (right). Solid, thin gray lines denote bounds for an energy 
consistent DWD given by assuming its final orbit must be Keplerian Dash-dot ted 
thin gray line marks the initial energy. For further details see lWoods et al.l (1201 ll) . 



-) 2 = (1 - y^r- 31 ) Vta + #)/(! + ( 4 ) 

where we use the notation qf = m^Jmd and q\ = m a /md. Taking af/a\ as a single 
parameter we can find the derivative of our solution with respect to the value of y: 



dcif/a; i 1 1 — qr 

J1L_L = _ 2B (A - yB), A = Vfei + <?f)/(l + <?i)- B = - —A (5) 

dy q f 1 + q\ 

where we have bundled up all of the mass terms in order to examine the effect of varying 
y, while holding other parameters fixed. For values (m^ = 2.0M Q , = 0.4M o , 
m a = 1.5M ) typical of an initial phase of mass loss in the formation of a DWD as 



prescribed by lNelemans et a l. (2000); lNelemans & Tou t (2005), we see that in the range 



1.5 < y < 1.75 we have -3.9 < d(af/ai)/dy < -2.5. Here a change in y by ~ 0.25 
leads to daf ~ da\. 

Clearly the fin al solution is q uite sensitively dependent on the value of y chosen, as 
was pointed out by Webbinki ( 20081) . and in fact an only slightly greater range allows for 



anything from merger events to greatly expanded orbits. It is then extremely important 
that the "correct" value of y be used for any given system. In figure 1, we see in the 
first panel that for a 2 + 0.5M o binary (with a 0.38M o hydrogen-exhausted core) y ~ 
1.505 would give a energy-and-angular-momentum consistent DWD system, whereas 
y* 1.75 would lead to a merger. If one were evolving this system forward, for example 
in a population synthesis study, the question becomes which is the appropriate value for 
the given masses. Similar behaviour can be seen for a system on the verge of an initial 
episode of mass loss (second panel). 
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While iNelemans & Tout! (120051) find 1 .5 < y < 1 .75 can provide plausible re- 
constructed initial conditions for all DWD systems, individual binaries can be fit with 
values of y ranging from 0.5 — 4 for given initial masses, differing from system to 
system. This relates to the fact that the final separation is very strongly dependent on 
the envelope mass available prior to the event. However, we have no means of know- 
ing a priori what the true initial masses, and so the correct value of y, would be for 
any known DWD we attempt to reconstruct. The difficulty remains that we have four 
unknowns: namely the initial separation; the separation after the first mass-loss event 
which stripped the envelope of the intially more massive component; and the initial 
masses of both components, while Nature provides us with only 3 constraints: the fi- 
nal separation, and the final remnant masses. This of course excludes any complicating 
factors such as the partial accretion by the older degenerate component of wind from its 
companion; any orbital evolution outside of the two main phases of mass loss (e.g. wind 
loss); or attempting to reconcile our understanding of the reconst ructed evolutionary 



pathway with the observed difference in ages of the components (ivan der Sluys et al. 



120061) . an uneasy task given our great uncertainty in determining the relative ages of 
white dwarfs. Without a physical theory underlying the y formalism which assigns a 
specific value of y for any given set of initial conditions, we cannot invoke angular 
momentum balance alone as a constraint upon the evolution of an arbitrary DWD. 

A lingering issue remains the relevant timescale associated with each episode of 
dramatic mass loss. While the second event must necessarily be a CE phase in or- 
der to produce the observed close orbit, the prior evolution is much more uncertain. 
In particular, the initial phase appears to eject the envelope with only a minor net 
reduction in the orbit, or even substantial expansion. A lthough typi cally presented 
as a d y namical timescale event, Ivan der Sluys et all (l2006h as well as iNelemans et all 
(2000): INelemans & Tout! ( 2005b only restrict themselves to timescales short compared 
to the nuclear evolution of the donor in their formulation, and in general a "y-prescribed 
event" could also be thought of as being a much longer, thermal timescale evolution- 
ary phase (Gijs Nelemans, private communication). However, an underlying physical 
mechanism for some spontaneous ejection of the envelope with little effect on the orbit 
remains lacking. 

For those DWD binaries in which we see evidence for an initial phase of orbital 
expansion with mass loss (namely, the older component being less massive), one would 
expect stable mass transfer (MT) to be a likely candidate. Indeed, though prior studies 
had claimed to have ruled this out, they did so under the assumption that any potential 
donor star which had evolved past the base of the giant branch would have too deep a 
convective envelope to allow for stable mass transfer, leaving a CE the only recourse 
upon reaching Roche lobe overfl ow (RLOF). However, we not e that MT is increasingly 
stable with greater core masses dHjellming & Webbinklll987h and that, particularly if 
one allows for non-conservative MT, stability is in fact quite possible. 



2. Stable MT followed by a CE 
2.1. Stability Criteria 

In order for mass transfer to proceed in a stable manner, we require that the response of 
the donor star to mass loss be such that it remain within its Roche lobe. We can define 
a response exponent £ such that R oc M^, which gives us the requirement £rl < £ a d 
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for stability against dyn amical timescale mass loss (where ^ is the adiabatic response 
of the donor, see iHjellmina & Webbink 1987) and £rl ^ <Tth> where £ t h is the thermal 
response of the donor (iSoberman et all 1 19971) . £rl is of course the response of the 
Roche lobe to mass loss from the donor; one can easily derive: 



c?log/?RL d In a dl^R^/a) dlnq 



where 7?rl is the Roche lobe radius. Through £rl's dependence on the mass ratio, 
the stability of mass transfer from the donor depends on the conservation factor /?, the 
fraction of mass lost by the donor which actually lands on the accretor. During case 
B mass transfer a long slow nuc lear timescale pha se is preceeded by a rapid mass loss 
phase on the thermal timescale ({Shore etlii1ll994l) . For our study we will take f$ = 0.3 
assuming that significant mass is lost from the system during this first stage of mass 
loss. The only non-compact objects we see to accrete at anything comparable to the 
thermal timescale are Herbig-Haro objects, in which we see the formation of bipolar 
outflows, and so we assume this is the mechanism by which matter is lost from the 
system here, as the accretor remains always well withi n its Roche lobe t hroughout 
stable mass transfer for the binaries considered here (see[Woods e t~al]|20lTI) . Mass is 



then lost carrying the specific angular momentum of the accretor. We find that this 
generally allows for stable MT for 0.83 < q < 0.92 and primary masses 1.0 - 1.3M Q . 

2.2. Double White Dwarf Evolution 

We can then evolve a set of 1.2+1. 1M binaries with varying initial periods through an 
initial phase of stable non-conservative case B mass transfer using a version of Eggle- 
ton's d etailed ste ll ar-evo lution code ev (also referred to as STARS, Eggleton] ( 1971 . 



19721) : iPols et all (I1995IV ). This leaves us with a set of systems in which the initial 
primary is stripped of its envelope, leaving a proto-helium WD remnant, an expanded 
orbit, and a companion which will necessarily enter into a CE phase upon its reach- 
ing the giant branch and RLOF, due to the dramatically upset mass ratio. We can then 
determine the result of a CE event by computing the Roche radius and the appropriate 
core mass for the companion to have reached that size, as well as using eq. 1 with 
®ce = 0.5. Doing so we find we are able to produce a string of models well in line 
with the observed distribution of DWDs whose older companion is the less massive, 
and in which neith er component has reached the helium flash (see fig. 2, as well as 



Woods etal. 2011) 



While promising, stable non-conservative MT as followed by a CE can at least 
explain a subset of the observed population of low-mass DWDs. The extreme effec- 
tiveness with which a CE phase removes angular momentum from the system means 
evolution through two such instances appears unlikely to allow for any binary to avoid 
merger without very careful fine tuning. Future studies will need to address other po- 
tential channels for DWD formation, in order to account for the entire population. 
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Figure 2. Period and component mass for four example binaries evolved from 
1.2+1.1M© systems with varying period, plotted just prior to the onset of a CE (pe- 
riod and first-born WD mass, circles) and after the CE phase (period and mass of 
younger WD, end of arrows). Known periods and masses of younger WDs are 
shown for single-line d (tr iangles) and double- lined systems (asterisks) as summa- 
rized in lMarshl d2011l) and lParsons et alj d201 lb . 
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